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The anisotropic Fermi superfluid of ultra-cold Fermi atoms under the p-wave Feshbach resonance 
is studied theoretically. The pairing symmetry of the ground state is determined by the strength of 
the atom-atom magnetic dipole interaction. It is k z for a strong dipole interaction; while it becomes 
k z — i/3k v , up to a rotation about z, for a weak one (Here f3 < 1 is a numerical coefficient). By 
changing the external magnetic field or the atomic gas density, a phase transition between these two 
states can be driven. We discuss how the pairing symmetry of the ground state can be determined 
in the time-of- flight experiments. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.90.+q 
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I. INTRODUCTION 

The trapped Bose or Fermi atomic gases pj are 
usually weakly interacting since the inter-particle dis- 
tances between the atoms are typically much larger than 
their scattering lengths. However, it has been recognized 
that the inter-particle interaction can in fact be tuned 
via Feshbach resonances 0, Q . The interaction changes 
from weakly to strongly attractive by varying the exter- 
nal magnetic field across the Feshbach resonance. For a 
two-component Fermi gas, the ground state evolves from 
a Bardeen-Cooper-Schrieffer (BCS) superfluid with long- 
ranged (compared with inter-particle distances) Cooper 
pairing to a Bose-Einstein condensate (BEC) of tightly 
bound pairs d, 0, • This cross-over has recently been 
a subject of intensive theoretical @, 0, 0, 0, 0] and 
experimental 0,0,0,0,0 investigations. 

Most previous investigations deal with s-wave Fesh- 
bach resonance. Therefore, both the Cooper pairs and 
the bound states between two fermions have s-wave sym- 
metries. Recent experiments demonstrated p-wave Fesh- 
bach resonances Qui inim. It thus raises the pos- 
sibility of p-wave Fermi superfluid and BEC of p-wave 
"molecules" jHHEl ( see also E3)- 

To begin, we first recall the well-known superfluid 3 He 
|2ri|. 3 He has (nuclear) spin 1/2, in general not spin- 
polarized and basically spatially isotropic, that is, the 
interaction between two atoms is dependent of the ori- 
entation of their relative distance. Back in the 60's, An- 
derson and Morel (2?| investigated theoretically the su- 
perfluid state for this system by assuming that the pair- 
ing exists only between the same (say f) species. They 
showed that the ground state corresponds to Cooper pair- 
ing in the I = 1, m = 1 channel, that is, the pair- 
ing wavefunction has the symmetry of the spherical har- 
monic Yii(fc) oc [k x + iky) (or its spatial rotation, such 
as k z — iky). This orbital symmetry is realized in the 3 He 
A-phase and is known as the "axial" state (2|| . One can 
thus suspect that for a spin-polarized but otherwise spa- 
tially isotropic system, the pairing is again in the Yn(fc) 



state. This has been confirmed in |2J|. 

However, in the atomic gas system, as demonstrated 
and explained by Ticknor et al. |l9| in the context of 
40 K, the magnetic dipole interaction, which breaks rota- 
tional symmetry, may be important. (An Alkali atom has 
one additional electron beyond closed electronic shell(s), 
hence it possesses a magnetic moment mainly due to this 
extra electron, though the precise value of this moment 
depends on the hyperfine interaction and the external 
magnetic field). They show that, due to this magnetic 
dipole interaction, for magnetic field along z, the 1 = 1, 
7tl = resonance occurs at a hig her magnetic field than 
the m = ±1 ones (see also |21|) For a given magnetic 
field, the induced effective interaction is anisotropic. In 
our previous paper (see also HH), we discussed the- 
oretically the expected ground state order parameter un- 
der this circumstance. Since the interaction of the m = 
channel is more attractive than that of the m = ±1, 
in the case that the m = and m = ±1 resonances 
are sufficiently far apart, we showed that the pairing can 
only occur in the m = channel, and the ground state 
becomes Yi (k) oc k z . The system is either the BCS 
state with k z Cooper pairing or the BEC state of bosonic 
"molecules" in the m — channel, depending on the de- 
tuning. Fcrmion pairing or bosons in m = ±1 channel 
do not exist. The orbital symmetry of the pairing is the 
same as the "polar" phase in the 3 He literature |2fij . 

For the case that two resonances in m = and m = ±1 
channels are closed to each other, we showed that the 
ground state symmetry is k z — i(3k y up to rotation about 
z, where z is along the magnetic field direction and f3 < 1. 
This state is thus intermediate between the axial ((3 = 1) 
the the polar (/3 = 0) phases. 

This pap er provides some details of our earlier investi- 
gation 28], as well as some additional information. It is 
organized as follows. We first study the Cooper pairing 
symmetry of a weak-coupling BCS model, where the p- 
wave pairing interaction is anisotropic. We show in that 
case how the analogous results cited above arise. Then 
we return to the slightly more involved case of Feshbach 
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resonance and determine the general phase diagram for 
this case. We first provide the results for the case where 
the system can be described by an effective single channel 
model, i.e., purely in terms of atoms interacting with each 
other via a two-body instantaneous interaction. Then we 
study the necessary modifications to these results when 
this approximation is relaxed. Finally a probe of the 
pairing symmetry based on the time-of-flight experiment 
is discussed. 



II. LESSONS FROM ANISOTROPIC BCS MODEL 

In the spirit of the work of Anderson and Morel 
|27). we first consider a weak-coupling BCS model of a 
Fermi system with one spin species under the interaction 
—Vt_%i ■ (We shall however study the case where this V is 
anisotropic, see below). The Hamiltonian of our system 
is H = Hf + Hy, where 



H 



f 



(1) 



k,k' ,q 



(2) 



Here e k = h 2 k 2 /2M, and fi is the chemical potential, and 
is the annihilation operator for a Fermion with wave- 
vector k. Standard BCS theory proposes the following 
trial ground state wavefunction 



iG)=nK-+^4 a - fe -)i°) 



(3) 



with normalization condition \u^\ 2 + \v^\ 2 = 1. Minimiz- 
ing the ground state energy (G\H\G) with respect to iig 
and i>£ gives the excitation spectrum 



where the energy gap satisfies 

^k ^ ^ ^k — k' 



2Ep 



(4) 



(5) 



In case the interaction is p-wave, only the I = 1 compo- 
nent survive in Vg_g,. We write 

m— 1 

V %~k>=^ E V rn {k,k')Y lm {k)Y* m {k') (6) 

m=— 1 

We shall study the situation where V\ = V-i < Vq. We 
also decompose the energy gap in spherical harmonics of 



I = 1 partial waves, 

Ag = A/(fc) 

TCI— 1 

= A frrrXlmik) 



(7) 
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with normalization J2 rn \fm\ 2 — 1- Substituting Eqs.®- 
into Eq.©, we have 



/(fc) = 7r^^y m (fc,A;')F lm (fc) 



[(£*,- A1 ) 2 +A2|/(fc')| 2 ]^ 



(8) 



Suppose further that V m (k, k') = V m is independent of k 
and k' , Eq.JSJ becomes 



f m = irVm E 



We perform the integration in fc in the following way, 

"A 



(9) 



E 



Po 



A 



de / dfl 



(10) 



where po is the density of state around the Fermi surface, 
A is the energy cutoff, and dfl is solid angle. Then we 
get totally 3 nonlinear equations, 



2A 

/ m io g (— )- / dnY* m (k)f(k)iog\f(k)\ 



(11) 



for m = — 1, 0, 1. 

We would like to find the solution to Eq. Ijllfl with the 
anisotropy characterized by the ratio Vi/Vq. Particular 
solutions to Eq. ljTJ can be easily found at two limits. 
When Vi/Vq = 0, oc Y 10 (k) oc k z because no other 
pairing can be possible. At the other limit, Vi/Vq = 1, 
the system is again isotropic the ground state pairing is 
Ajj oc Yn(fc), or its rotations. In terms of the language 
by f m % for Vi/Vo = 0, /o = l,/i = f-x = (up to 
a gauge transformation). For Vi/Vq — 1, a particular 
solution is /o = l/\/2, f± = f-i = 1/2 such that 



v/2 



(12) 



Note that it does not contradict with the pairing Yn (k) oc 
k x + iky since these states are degenerate in an isotropic 
limit V\ = V-i = Vq. With gradually tuning the 
parameter < Vi/Vq < 1, we will show below that 
A(fc) oc k z — i[3k y up to a z-axis rotation with j3 < 1. 

Before solving the gap equation Eq. (|ll|l for the general 
intermediate values of the pairing < Vi/Vq < 1, it is 
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helpful to first identify the critical V* /Vq such that /i 
and /_i start to deviate from zero (i.e., when the pairing 
starts to deviate from Y w )- Let us take the gauge where 
fo is real. Linearizing Ea. l|ll(l in fx and /_i around the 
point (f , f u = (1, 0, 0), we have 



_W h 

2tt P oV 1 * \ /*! 



Ai A 2 
A 2 A x 



h 



(13) 



where A x = log(2A/A) - J dn\Y n \ 2 log |y 10 | - 1/2, A 2 = 
1/2. Thus when we increase Vi, the first non-trivial so- 
lution for fx and /_i is obtained when (fx, fix) gives 
the largest eigenvalue of the matrix. It can be easily 
seen that this corresponds to the solution fx = fix, with 
eigenvalue Ax + A 2 . The critical value V* thus satisfies 



2A f 

io g (— )- y dn|r n | a iog|y u 



, ,, _ . „„ (14) 

2tt Po V 1 * |_ A 7 J 

On the other hand, at this point (from Eq. I|ll() with 



in 



0), 



OA f 

io g (— )- y df>|y 10 | 2 iog|y 10 | 

Eliminating log(2A/A) in Eqs.{TIJ|-||T5J| gives 



2irp V 



Vx*/V 



1 



1 + 2np Vo 



(15) 



(16) 



For V* /Vo < Vx/Vo < 1, f±x becomes finite. Provided 
no additional phase transitions occur, we expect then 
the relation fx = /*j remains valid. Furthermore, our 
system is rotationally invariant around z. Under this 
rotation, fx — > f\e % ^ but /_i — > /-ie - *^ where is the 
rotational angle. By suitable choice of <j> we can thus 
make both f x and /_i = /* real and positive. We shall 
confine ourselves to this case without loss of generality. 
Under this choice the order parameter has the form k z — 

Hence, when V* /Vo < Vx/Vo < 1, we search for the 
solution fo > 0, fx = f-i > under the normaliza- 
tion condition f§ + 2ff — 1. The solution fo (where 
fx = f-x = v/(l- / 2 )/2 ) for different Vx/V and p V 
is presented in Fig^ For given poV), the pairing is k z 
(fo = 0) at small Vx/Vo until the latter exceeds the criti- 
cal Vx/Vo- Beyond this critical Vx* /Vo, fo decreases grad- 
ually with increasing Vx/Vo- At Vx/Vo = 1, fo = 1/V2 
for all poUo's meaning that the pairing becomes k z — ik y 
(corresponding to f3 = 1). 



III. P-WAVE FESHBACH PAIRING 

We return to the case of p-wave Feshbach resonance. 
We begin with the total Hamiltonian H = Hf +Hb + H a 
where again Hf is for the free fermions, 



H 



f 



k 



(e k - n)ata % . 



(17) 




FIG. 1: 

legend. 



fo as a function of Vx/Vo for PqVq's given in the 



H\, is for the molecules in the "closed" channels, with b^ m 
the annihilation operator for a Boson with angular mo- 
mentum I = 1 and z-axis projection m with momentum 

q, 



(18) 



Here (5 m 's are the (bare) detuning of the I = 1, m res- 
onance. H a represents the Feshbach coupling, which is 
a gene ralization of the ones already commonly employed 
[8J, |30j for s-wave Feshbach resonances to the case of sev- 
eral I = 1, m = 0,±1 closed channels ESI El: 



Vm (^) b {m a -k + q/2 a k + q/2 



h.c. 



(19) 



where <p m (k) = — ^gyf kYx m (k)a m . The factor Yx m (k) re- 
flects the symmetry of the I, m bound state and the linear 
factor in k arises from the small momentum approxima- 
tion for the coupling. a m is the corresponding coupling 
constant, it has dimension [i?][X] 5 / 2 . More precisely a m 
can be an analytic function of k. For low energy (den- 
sity) phenomena, we shall however only need its value 
a m = oi m (k = 0). 

We shall apply mean-field theory to the Hamiltonian 
H . In this approximation, we discard all terms that in- 
volve finite momentum Bosons b^ m with q ^ 0, and re- 
gard b^ = o. m as c-numbers. It is then convenient to intro- 
duce 



D„ 



a m bo. 



(20) 
(21) 



A£ = ^D m kY lm (k) 

771 

Hence the mean-field Hamiltonian for the Fermions be- 



comes 



H°* = - n)a\a % + (A^o_gog + h.c.) (22) 



k 



This can be solved via the familiar Bogoliubov transfor- 
mation, which gives 



(23) 



Minimization with respect to the c-numbers b^ =0 nl give 



b 0>m (8 m -2n) = -Y] ip* m {k){a_zat) 



k 



which can be written as, 



{5 m ~ 2fj)D Tl 



2tt- 



\a m \ 2 kY* m (k)A, 



L3 I [(e k -^ + \A^\ 



(24) 



(25) 



The R.H.S. of this equation is formally divergent (This 
divergence is only formal: the sum may actually converge 
if the k dependences of a m (k) are included) Nevertheless, 
it can be " renormalized" by re-expressing this equation in 
terms of parameters that enter the two-body scattering 
amplitude (c.f. |2J[and also the theoretical references in 
the s-wave case |Ml^ llollllTll2l |.) Let us then consider the 
scattering of two particles in vacuum. We denote ftik') 
as the scattering amplitude for two incident particles with 
relative momentum k and out-going momentum k'. We 
can define f m (k) (do not confuse this f m {k) with those 
in Sec II) via 



/*(*') = E f m {k)i^)Y lm {k')Y* m {k) 



(26) 



This equation is the generalization of the corresponding 
one in scattering theory to the case of anisotropic inter- 
action (but still with rotational invariance around z). It 
reduces to those in standard quantum mechanics text- 
books |3l| if fm(k) is independent of m. fs{k') is related 
to the T-matrix by 

ML 3 , 



AttK 2 



WS = 2e fc +iO+) 



(27) 



The T-matrix can be evaluated in the standard way: 

rn 



E 



<lm( fc ')lm(fc) 

E-8 m -U.W(E) 




200 

B (G) 



FIG. 2: Plot of the fit to the experimental results of 
Ref . for —1/vm and — c m for 40 K atoms in the |/, m/) = 
j 9/2 , —7/2) hyperfine state, vi = w_i and ci = c_i. — l/u m 's 
vanish at the resonant fields _B„, with Bq > B±i- 



Q f (p, uj) is the Greens function for free Fermions at 
wavevector p (p has dimension of [L] -1 ), frequency u>, 

and Of,™** (E) i s the Greens functions of free TO-Boson at 
zero momentum q = and energy E. Here we have al- 
ready made use of a simplification due to the fact that 
Tj( mm ) vanishes for m ^ m' (see Eq. (|2*9"|) 1. hence there 
are no cross-terms between different m's in Eq. 128|) . 

Substituting Eqs.®-© in to Ea.l(27)). and using the 
definition l|26|) . we get 



M 



2e k 



4ttS 2 f m (k) 



r.3 



L 3 ^ 2e k - 2e„ + ^0^ 

p 



(30) 

The sum over p in the R.H.S. formally diverges. However, 
we are interested only on its dependence on k at small 
k's. We thus expand the R.H.S. in powers of k: 



2e k 



1 



P 



L 3 2e k - 2e p + 



a, 



■__MJ_y 

I 2 h 2 L 3 ^ ' 



+( 



Ma r 



M 1 



h 2 l 3 ^ p 

v 



Ei) fc2 



M 



fc 3 (31) 



where 



In the low-energy limit, Ea. H3U|) - (|31|l should reduce to 
(28ihe general result for the low-energy scattering between 
a pair of particles. Using the same notation adopted in 
H.ef. [lflj ■ the scattering amplitude is parametrized as 



p,LJ 



fm(k) 



Cmfc 



ik 3 



(32) 



E 



ym(P>m'(p) 

E — 2e„ 



(29) 



where the dimensions of v m and c m are [L] 3 and [L] 1 , re- 
spectively The magnetic field dependent parameters v m , 
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c m are in principle available experimentally. An example 
is as shown in Fig[2J Our renormalization scheme is to 
absorb the bare parameters S m and a m into the physi- 
cal (renormalized) parameters v m , c rn by identifying the 
equivalence between Eas. (|3U|) - (l31ll and Ea. l|32|) . Hence 



1 



4irh 2 8„ 



--YD 



Anti 2 



M M\a r , 



M 1 



P 



(33) 
(34) 



In Eq. l|3"3")) and Ipjljl. if the p dependence of a m 's 
had been included, the sums involving p should have 
the extra factors |a m (p)| 2 /|a m (0)| 2 . When treating the 
many-body problem, as we shall see below, the only 
"divergent" sums are again |a m (p)| 2 /|a TO (0)| 2 and 
Eff\a m (p)\ 2 /\a m (0)\ 2 p 2 . We shall apply Eq. ^ and 
()34JI to eliminate these sums in favor of v m and c m . Af- 
ter this procedure, low energy phenomena (involving only 
small p's) can be parameterized entirely by v m , c m and 
a m (0). To simplify our presentation, we shall not include 
explicitly the factors |a m (p)| 2 /|a m (0)| 2 in our calcula- 
tions below. 

Returning to Ea. (|25[) . we rewrite its R.H.S as 



where 



h(k) = 



[( £fc - M ) 2 + |A fc H 2 ] 



1/2 



k 2 , fl 
+ — 



|A fc -|\ 
2e 2 > 



(36) 



2?r x - 



2 kY* m (k)A, 



I [(e k -ri 2 + \AtfY 



2TT\a m \ 2 y D m ,Y? m {k)Y lm ,(k)k 2 



Y,D m .Yt m {k)Y lm .{k) 



L 3 



xU(jfe) + — (1+ M 



|A fc H 2 , 
e k 2d > 



(35) 



The sum involving h(k) in Eq. 1|35[) is convergent. For 
the rest of the terms, their divergences are the same as 
in either Eqs. (|33|l and (|34fl . These sums can thus be 
expressed in terms of the physical parameters v m and c m . 
Using Eq. (|35|l , the " gap equation" Eq. (|25[) can therefore 
be rewritten as 



= L3~Y, D m Y* {k)Y lm {k)h{k) (37) 



= |5^£» ro y 1 * 1 (^imW/i(fc) (38) 



L 3 , 

A: , m 



and a corresponding equation with to = 1 <-> m — 



-1. 



Ea. (|37() - (|38() are to be solved under the constraint from 



6 



the number equation 



(39) 



The first term is from open-channel atoms whereas the 
second term is from the closed-channel molecules. From 
BCS theory, 



1 



[( efc - M )a + |AgP]i 



The first term in Eq. 1(39(1 is thus again formally diver- 

V**"' 



gent: (atar) <~ at large fc. However, this sum it can 



be treated in analogous manner as Eq. 1)35(1 by employing 
again Eq. The result is 



-to 



1 



[(e* ~ M) 2 



47T ^ 47T?l 4 



M 2 C. 



|A fc H 2 ] 



I -Dr. 



24 } 



-T 



I A, 



(40) 



Eos.lpTTjl. J2H1), and ijjQjl are our principal equations, 
with parameters characterizing the Feshbach resonances 
entirely in terms of v m , c m and l/|a m | 2 . These equations 
determine the order parameters D m and chemical poten- 
tial /X for given density n and "interaction parameters" 
v m , c m and l/\a m \ 2 . 



IV. BEC-BCS CROSSOVER AND QUANTUM 
PHASE TRANSITIONS 



into account higher order terms in k 2 in denominator of 

fm(k). 

Moreover, as explained in Ref.fljj, due to the dipole 
interaction, B'q > B\ = B*_ 1} as can be seen again in Fig 
121 Thus in the field range of interest, — l/v\ = —1/v-x > 
— 1/vo. We can say that, at a given field, the effective 
interaction between the Fermions is less attractive for 
relative angular momentum projections m = ±1 than 
to = 0. 

Since the interaction is less attractive for angular mo- 
mentum projections to = ±1, for sufficiently large differ- 
ence between —1/vq and —l/v±i we expect (and verify 
below) that the pairing is entirely in the to = partial 
wave. We thus first begin our analysis by assuming that 
only Dq is non-vanishing. 




Generally with Feshbach resonance for the sub-channel 
to, l/v m is field dependent, vanishing at the resonant 
field B m . In contrast, c m has a definite sign (c.f. Fig. 
|2J). For the ease of discussions, we shall assume that 
Cm < and field independent, —l/v m is an increasing 
function of field —l/v m > (<)0 for B > (<)B m , as in 
the case of 40 K. (This corresponds to the case where S m 
is an increasing function of field and a m weakly field 
dependent, c.f. Eas.ffify and (J2}). For B < B m , a 
bound state appears. The energy of this bound state 
is given by — ej, im = ~h 2 K m 2 /M with k = in m being a 
pole for f m (k), that is, —l/v m - Cm^m - n m = 0. For 
given to and small detuning below that resonance, n m is 
small and is given by K 2 n — l/[(—c m )(v m )}. Since l/v m 
should be roughly linear in B near the resonance, eb iITl 
increases linearly with (B m — B) (in contrast to s-wave, 
where it is quadratic; see also |24j) These results for small 
detuning apply provided K m « (— c m ), or equivalently 
l/w m <C (—c m ) 3 . We shall always confine ourselves to 
this regime for negative detunings (even when we write 
B <C B m ). For larger negative detunings we need to take 



FIG. 3: The dimensionless parameters Do and jl as functions 
of —1/vo- £o = 5i = —100 and D±i = in this case. The 
dashed line represents —e b /2tF- 

For simplicity, we shall first drop the terms with ex- 
plicit l/|a m | 2 factors. If this is valid, then our equations 
(|3Tjl . P5)l. (|4*U|l are same as a corresponding model con- 
sisting only of Fermions interacting among each other 
with an interaction such that the corresponding two- 
body scattering amplitude parameters are given by the 
same c m and v m . In analogy with the s-wave literature, 
we shall call this the single-channel approximation (note 
there are still three (to = —1, 0, 1) Feshbach resonances). 
The effect from finite l/\a m \ 2 will be discussed at the 
end of this Section. 

With these simplifications, Eq. (|37|l and 14(J|I become 

(41) 
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and 



A- 

M 2 



[(6 fc - M )2 + |A fe -P] 



2e 2 j 



(4tt) 2 /i 4 



(-co)l^ol 



(42) 



respectively. Eqs.JSJJ and (|40|l can be solved simultane- 
ously similar to the s-wave case. By gauge invariance, 
we shall choose Do to be real without loss of general- 
ity. It is convenient to express the results in dimension- 
less form. We thus define jl = fJ,/ep, D m = D m /hvY, 
c m = 7i _1 / 3 c m , and v m = nv m where ep = h 2 k F /2A / I, 
vp = hkp/M, and k F = 6ir 2 n. The results are as shown 
in Fig. |3](for the case Co = c\ = —100, see below for the 
reason of this choice) . 

In the BCS regime, corresponding to "large" external 
magnetic field B ^> Bq thus — Vq 1 — > oo, we have fl ~ 1, 
and Dq ~ 0. We can ignore the term oc \Dq\ 2 on the 
left hand side of Eq. (|41(l . This equation then reduces to 
the corresponding one in Sec II with the effective pairing 
interaction Vq in the m — channel there proportional 
to (-l/« + 2Mc of i/h 2 )- 1 = (-1/vo + cofc 2 ,)" 1 . This 
combination follows from the factor (So — 2/z) on the left 
hand side of Eq. lj2~5|) . 

In the BEC regime, corresponding to B <C Bq or 
—Vq 1 — > — oo, /i ~ — £6.o/2, and Do approaches to a con- 
stant. In Eq. I|41|) we can ignore all terms explicit in A^ 
or Dq. Writing k = (2A/|^|) 1 / 2 and performing the sum 
over fc, we find that k obeys —1/vq — cqk 2 = k 3 . Com- 
paring this with the equation for the two body bound 
state, we find = e;,o/2 as claimed. In Eq. $£2$ , 



we can expand the square-root [(e^ — fl) + |Ap| ]" 
(e k + - |Ar| 2 /2(e fc + \fi\). We then get 



(4tt) 2 H 4 



\Dq 



(-co) - 2 K 



-1/2 



(43) 



where k was defined above. For small detuning k <C 
(-co), i.e. l/do <C | Co | 3 , a relation well-satisfied in the 
range in Fig |21 and others below, we then obtain Dq — 
(32tt/3) 1 /3(- 5o )-i/2. 

The "cross-over" behavior in Fig|3] is analogous to 
the s-wave case, where the corresponding x-axis is x = 
— l/(n 1 / 3 a s ) where a is the s-wave scattering length. 
Note here Do has the dimension of [E][L] _1 and behaves 
differently from the s-wave A in the BEC limit: it satu- 
rates rather than continue to increase. Rigorously speak- 
ing, we actually expect a (quantum) phase transition at 
the point where \i crosses zero 0, E3 : the system is gap- 
less for /i > but gapful if /i < 0. However, we found no 
evident changes of slope in /i and A when this transition 
is crossed, and thus expect this phase transition is very 
weak thcrmodynamically. Similar results were found in 
Ref H3| for the axial (D x ^ 0, Dq = D_ 1 = state). 



We have also performed calculations for other values 
of Co- The size of the crossover region is roughly pro- 
portional to the value of cq (which is in turn inversely 
proportional to the density). For example, for Co = —200 
(not shown), the corresponding results can be captured 
well by replacing the x-axis by — l/2vo and dividing Dq 
by l/v^in FigH 

The above behavior applies only to sufficiently large 
— — (—Vq 1 ) > 0. When this difference is sufficiently 
small, D±i will become finite. Similar to the treatment in 
Section II for the anisotropic weak-coupling BCS model, 
we first identify the critical v^ 1 such that D\ and/or 
D_i first start to deviate from zero. For simplicity, in the 
results below we shall take Co = c±\. This is satisfied to a 
good approximation at least for 40 K near the resonances 
(see Fig EJl . We linearize Eqs.(|2Z|l-(|SHl) (and the latter 
with m — 1 <-> m = —1) in D\ and D_\. We obtain 



M ( Dx 



Airvl \D*_ 1 J \A 2 Ax J \ D*_ x 



Ax A 2 \ ( D 



(44) 



where 



A > = ^J2\Yn\ 2 h(k)-^J2 D o\ Y ^\ 2 \^\ 2 9(k) 

k k 



M 2 cx, 3M , 



2ttH 2 ^ lOirfr 2 

A2 = jE^oHttiiy*) 

k 

k 4 1 



(45) 



1 



2 [{e k -ri 2 + D 2 \Yxo\ 2 k 2 f 2 4 



) (47) 



Here h(k) is as defined in Eq. i|3(3|) but with only Dq ^ 
0. We can verify that A 2 > 0. Therefore the largest 
eigenvalue for the A-matrix, hence the smallest value for 
v*, belongs to the class Dx = D*_ l . The corresponding 
critical value is then given by -j^jv =Ax + A 2 . On the 
other hand, from Ea. l|37(l . we have 



M 



2n 



XQ\ 2 h(k) 



M 2 



Co 



27r^ 2 



(A*- 



-^Ld 2 ) 

20nh 2 o) 



(48) 



Combining these two equations and re-writing it in di- 
mensionless form, we have 



1 



3(6tt 2 )! 



D 2 o(-5 ) 



v* vo 5ir 

n f°° , f 1 , x 4 (l-3y 2 ) 
+9tt / dx / dy i — t-t - 

Jo J -i [(x 2 - fi) 2 + ^D 2 x 2 y 2 Y/ 2 

(49) 

In the BCS limit, D 0, the first term in the R.H.S. 
of Eq. H49[) is negligible whereas the second, which we shall 



call Ki, is finite. The dominant contribution to K\ oc- 
curs near x ~ yffi, where the integrand diverges if Dq 
were zero exactly. To evaluate K\ in the limit Dq — > 0, 
we add and subtract the integral 



K 2 = 9tt 



r7i 



dy 
-i 2[£* 



(l-3y 2 ) 



2„21 



(50) 



where ^ = x 2 — jl, and £o 3> -Do is an arbitrary cut-off. 
The integrand for K\ — K 2 now has no divergence near 
x 2 — u sa 0, and their difference vanishes in the Dq — ► 
limit due to the y integral. X 2 can be evaluated by first 
integrating with respect to £, we get 




l/v„ 



-18tt y (l-3y 2 )logy 
12-k 



(51) 



In going from the first to the second line, we have used the 
fact that J^dy (1 — 3y 2 ) = 0, hence the terms involving 
£o and Dq drop out in the £q ^> Dq limit. Therefore in 



the BCS limit we obtain -5! 



12tt = 37.7. 



In the BEC limit, recalling that Dq approaches a con- 
stant ~ (327T/3) 1 ' 3 (— Co) -1 / 2 whereas u is large and 
negative, we find that the second term in the RHS of 
Eq. (|49|) is negligible compared with the first (provided 
k <C (-co)). Using the value for Do, Eq. f$ty yields 
— v*^ 1 + v^ 1 — ► 487r/5 = 30.2. In the intermediate 

as a func- 



-v. 



„~.-i 



v, 



regime, the numerical results for 

Vq 1 is shown as the thick black line (Dq = 0) in 



n 



tion of — v n 
FigH 

We now solve for the order parameters D m when 
— 1 + Vq 1 is less than the critical value. As already 
mentioned, we have already made use of gauge invari- 
ance to choose Dq to be real. Under this choice, the 
solutions we found belong to the class D\ = D*_ v Writ- 
ing D\ = \Di\e lx , then has the angular depen- 
dence oc DqIc z + i\/2\Di\k ■ a cx (k z — if3k ■ d) where 
a = (cosx)y + (sinx)i is a unit vector perpendicular 
to z and (3 = \/2\Di/Dq\. A particular solution is given 
by the case where D\ and D_i are both real where a = y. 
Other solutions are simply related to this one by a rota- 
tion about z. Without loss of generality, we shall there- 
fore assume that D m 's are all real below. 

The contour plot of (3 (by solving Eq. P7 |l -l|30 |) with- 
out the l/|a m | 2 terms) is shown in Fig0| (see Fig. 2 of 
|28| for a plot of D±i) (3 monotonically increases down- 
wards or towards the right in this diagram. In the (3 = 
(D±i = 0) phase, the state is rotationally invariant about 
z, whereas this symmetry is broken in the (3 ^ phase. 
There is a (quantum) phase transition between these two 
phases when one crosses the critical line — v*^ 1 + v^ 1 . 



FIG. 4: Contour plot of (3 as a function of — 1/vi + l/«o and 
— 1/vo for co = c±i = —100. The line for j3 = corresponds 
to the critical value — 1/v* + 1/vo- 




FIG. 5: Contour plot of (3 as a function of — 1/vi + l/«o and 
— 1/vo for Co = c±i = —200. 



For co = —200, the results are qualitatively the same 
provided we double the value of —X/vq. The correspond- 
ing plot is shown in Fig[SJ 

In the intermediate splitting regime, we thus predict 
the state is k z — i(3k y on the BCS side, whereas it is 
k z on the BEC side. For large positive detuning, the 
splitting should be less relevant and the pairing state 
should resemble more that of the isotropic system. On 
the BEC side, the system should be closer to a Bose 
condensate of lowest energy molecules (k z ). 
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For the 40 K case studied in Ref . fl9| . the Feshbach res- 
onances are at Bq ps 198. 8G and B{ ps 198.4G. There, 
C\ and Co are both only weakly field dependent and are 
approximately given by — 0.02a " 1 (see Fig|2J). Our choice 
of c = — 100 above corresponds to a density of roughly 
10" n ao 3 = 6.7 x 10 13 cm- 3 . The range of -1/vq for p, to 
go from to 1 corresponds to roughly 0.5G if we take the 
same gas density. (This field range is proportional to n 2 / 3 
). Near the resonant fields, — 1 + v^ 1 ps 2.1 x 10~ 8 aQ 3 
and is roughly field independent. Thus the density de- 
termines the values for both co,i and — v^ 1 + v^ 1 while 
varying the magnetic field corresponds roughly to mov- 
ing along a horizontal line on our phase diagram of Fig0] 
(with increasing field towards the right and the distance 
of the line from the horizontal axis proportional to n _1 ). 
For the density cited above, — #7/ 1 + v^ 1 ps 2000, hence 
we expect only the k z phase to be observed. To observe 
the phase transition, we need a 40 K gas of a much higher 
density, or another gas species with much smaller split- 
ting between the m — and m — ±1 resonances. 



In the above treatment, we have assumed that the 
terms with explicit \/\a m \ 2 factors in Eas. l|37[l - i|40[l can 
be dropped. The validity of this assumption has to be 
investigated by first-principle calculations of two-atom 
collisions. Following Ref |37|, we define the quantity 



_ M 2 |a m | 2 ^ 1 



£ 4 L 3 



~ P 2 



(52) 



r\ m quantifies the relative strength of coupling between 
the m closed-channel and the continuum. With this no- 
tation, Eq. (fMj) becomes 



2U 12 \L + ri m ) 



M 2 a 



(53) 



The single channel approximation used above corre- 
sponds to r] m — > oo. Wc now investigate the effect of 
finite r] m on the many-body problem. Expressing the 
Eqs.lO, and (gDJ) by T] m instead of l/\a m \ 2 , we get 



^J2 D ™ Y *o(k)Yi m (k)h(k) (54) 



L 3 

k 



+ 2^fi, _ [ftPlftl" + 2|D |* + 6,^) - DUDS] 

47r v\ n z 1 + rji 407r^/i 4 L J 



— ^D m Y^{k)Y lm (k)h{k) (55) 



k.m 



2i 3Z ^ [(e fe - M ) 2 + |Ad 2 ]5 



fe 
M 2 



\*S\\ 
2e 2 j 



(4tt) 2 /i 4 



Y- 

^ 1 



»7n 



I A, 



(56) 



We note here that the last term of Ea. (|56fl involves a 
factor (2 + rj m )/(l + r/ m ). This term actually arises from 
two contributions. The first one, involving 2/(1 + T} m ), 
is due to the closed-channel molecules (the last term in 
Eq.|@DJ|), and a second one, involving r] m /(l + rj m ) due 
to the open-channel atoms. 



To show the effect of finite 770, we plot the general 
results for Dq and jl in the state D±i = in FigHj] 



The value of rjo mainly affects the results on the BEC 
side. In fact, deep in the BCS limit the results are 
independent of ry , as expected since all atoms are ba- 
sically in the open-channels, jl is still basically linear 
in l/vo with the same slope, and thus the width of 
the "crossover" region remains to be determined by Co 
only (but not 770)- hi the BEC limit, when rjo decreases 
from 00, the contribution from bound state molecules 
becomes more important (see in particular Eq. i|56|) 
and the discussions below it). The chemical potential 
jl becomes closer to the two-body value for decreasing 
7]q as shown in Fig. E{b). Dq has the limiting value 
~ (32^/3) 1 / 3 (-c )- 1 / 2 [(l + 77o)/(2 + %)] 1/2 , and thus 
decreases with decreasing 770. 

The phase transition between the k z state and the k z — 
i/3k y state can be investigated as before. The value for 
—v^ 1 + Vq 1 is not affected by the value of ?y m 's in the 
BCS limit. In the BEC limit, if t] Q = r)±i = 77, we obtain 
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FIG. 6: Plots of (a) Do and (b) fi as functions of — 1/wo 
at £o = c±i = —100, D±i — for different 770 's. rjo = co 
corresponds that of Fig The dashed line in (b) represents 
to the two-body result. 



-v 1 ' + Vq 1 — ► (30.2)?7/ (2 + 77). This critical value thus 
decreases with decreasing 77. Its behavior is as shown in 
Fig [7| The 7/ dependence of this transition line has also 
been investigated in Ref |29j . 



ter a s in Eq.JTJJ satisfies h 4 / (\a s \ 2 M 2 ) < 1/n 1 / 3 . This 
case is often also referred to the "wide resonance" regime. 
Though for both s and p-wave resonances, the single- 
channel approximation would be valid for sufficiently 
strong coupling (large a) between the closed and open 
channels, for s-wave the condition for a single-channel ap- 
proximation is always satisfied when the gas is sufficiently 
dilute. In contrast, for p-wave the condition 77 ^> 1 is 
independent of the density of the gas. This difference 
can be understood by noting that the dimension of a s 
is [E][L] 3 / 2 , whereas that of a m here is [E][L] 5 / 2 . The 
effective range r s for s-wave Feshbach resonance [32, 
is proportional to h 4 / '(\a s \ 2 'M 2 ) , and the additional term 
oc J2(^/p 2 ) m Eq. ffi^ does not arise. 

For the s-wave case, the terms " single- channel approx- 
imation" and "wide-resonance" have been used inter- 
changably. However, we should note that for our p-wave 
case, the condition 77 3> 1 is not necessarily in conflict 
with the fact that the the peak in scattering cross-section 
discussed in Fig.l of Ref.^| is narrow. This p-wave 
scattering cross-section is proportional to |/ m | 2 , which 
is given by 



I/77 



(57) 
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l/v„ 



FIG. 7: The transition line 
and the k z — ij3k y 
-100. 



1 +v 1 between the k z phase 
phase for 77 's in the legend. £0 = c±i = 



We conclude with a discussion on the condition 77 3> 1 , 
or equivalently H 4 / (a m M) 2 <C \c m \, for the validity of the 
effective single-channel approximation. It is worth noting 
that the corresponding condition for the s-wave Feshbach 
resonance superfluidity is quite different. For this latter 
case, the sing le-channel approximation near resonanace 
is valid if |32l I331 13^. I35I l3q the corresponding parame- 



For B > B* n , where l/v m < 0, the cross-section 
thus has a peak at finite energy E = h 2 k 2 /M where 
k 2 = l/(\c m v m \). The denominator can be written as 
c 2 „(fc 2 - fc 2 ) 2 + k 6 . For B > B* m , |/ m | 2 as a function of k 2 
then peaks very near k 2 = k 2 and has a width approx- 
imately given by fc 3 /|c m | = l/(|c m | 5 / 2 (-7j m ) 3 / 2 ) < k 2 , 
and therefore a very narrow (in energy) resonance. This 
is entirely a result of the special form of / for finite angu- 
lar momentum. This narrowness of the peak is therefore 
not directly connected with the assumption on the mag- 
itude of 77 made here. 

In addition, the condition for 77 3> 1 here is also not 
necessarily in conflict with the assumption of a "narrow" 
resonance in the context discussed by Gurarie et al |29l| . 
They defined a parameter 7 oc cv^M 2 ?! 1 / 3 , and called the 
case 7 <C 1 to be a narrow resonance. Their condition 
is fulfilled always at sufficiently low densities. Thus the 
resonance can be narrow in their language whereas the 
single-channel approximation can still be valid if n 1 / 3 <C 

S 4 /(«mM) 2 « \c m \- 



V. EXPERIMENTAL PROBE OF PAIRING 
SYMMETRY 

It was suggested by Altman et al. [3£| that the mea- 
surement of the atom shot noise in the time-of-flight 
(TOF) absorption image can reveal the properties of the 
many-body states, including the symmetry of the pair- 
ing state. Later the proposal was carried out by Greiner 



11 



et al. [33 for 40 K atoms under an s-wave Feshbach res- 
onance. As discussed in previous section, the pairing 
state induced by p-wave Feshbach resonance is either k z 
or k z — i0k y where can vary under a change in the mag- 
netic field or the gas density. In this section, we propose 
how can be determined in a TOF experiment. 

In the TOF experiment, one can measure the quantity 

G(r,-r) = (n(r)n(-r)) t - (n(r)) t (n(-r))t 

K ( n k n -k) ~ ( n k)( n -k) ( 58 ) 

where k = mr/(ht). t is the time after the trapping 
potential and the interaction between atoms are suddenly 
turned off. Using the mean-field Hamiltonian in Ea.l)22|l. 
we have 



Q(f,—r) oc 



±bt{kl + p~kl) 



(59) 



where k = k/kp. The two phases k z and k z — i0(k ■ a) 
can be distinguished by the presence of anisotropy in the 
x — y plane. For example, consider the weighted average 
over <f> defined by 



S (a) = (cos (2(0 + a)) 0(f, -r)) 



(60) 



This average vanished in the k z phase. For 0^0, So (a) 
becomes finite and can in principle be evaluated. The 
presence of the k z + 2 k 2 term in the denominator com- 
plicates the analysis. There are nevertheless two possible 
simplifications. In the BEC limit, the terms involving 
D 2 . in the denominator are negligible compared with the 
first. In the BCS limit we can confine ourselves to data 
where r 3> Rf, equivalently k ^> fcp. In both cases the 
correlation function has the angular dependence 



G{r,-r) ex (K+FK) 



(61) 



Then, for a = (cosx)y + (sinx)i as before, we get 

So (a) cx ~ Dl[3 2 cos {2{ X - a)) . (62) 

Thus So (a) will reveal the finiteness of as well as the 
value of x hence the direction of a (up to a sign). The 
values of however are still difficult to determine in this 
way, both because the need to evaluate the proportional- 
ity coefficient in So as well as the fact that experiments so 
far detect a correlation in s-wave substantially less than 
the expected theoretical value. However, we can consider 
instead 



51 ee (cos 2 (f • z)Q{r,-r)) 

52 = (sin 2 (f • z)Q{f, —f)) 



These give 



Si 

So 



3 + 2 
2 + A0 2 



(63) 
(64) 



(65) 



and hence the value of can be determined via 
3 - 2(Si/S 2 ) 



Z 



4(Si/S 2 



1 



(66) 



For the polar state, = and Si = 3S2/2. For the axial 
state, = 1 with Si = 2S 2 /3. 



VI. CONCLUSION 

In conclusion, due to the magnetic dipolar interac- 
tion between two fermionic atoms, the superfluidity in 
trapped gases induced by a p-wave Feshbach resonance 
is expected to behave very differently from the more fa- 
miliar example of 3 He. The symmetry of the ground state 
can be tuned by varying either the density of the gas or 
the applied external magnetic field. The phase can be 
a polar state or a state intermediate between the polar 
and the axial phase. We also propose how the pairing 
symmetry can be detected from an experimental investi- 
gation. 
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